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Abstract 

A brief review is given of the Continuous Spontaneous Localization (CSL) model in which a classical field 
interacts with quantized particles to cause dynamical wavefunction collapse. One of the model's predictions 
is that particles "spontaneously" gain energy at a slow rate. When applied to the excitation of a nucleon in 
a Ge nucleus, it is shown how a limit on the relative collapse rates of neutron and proton could be obtained, 
and a rough estimate is made from data. When applied to the spontaneous excitation of Is electrons in Ge, 
by a more detailed analysis of more accurate data than previously given, an updated limit is obtained on 
the relative collapse rates of the electron and proton, suggesting that the coupling of the field to electrons 
. and nucleons is mass proportional. 
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1. Introduction 



It is appropriate to discuss comparison of experiment to a theory with fundamental pretensions in a 
volume dedicated to Dan Greenberger whose own theoretical work of a fundamental nature has seldom been 
I far from testability. 

O^. In standard quantum theory (SQT) the statevector evolves in two ways. One evolution precedes 

smoothly via Schrodinger's equation. The other is the abrupt (and ill defined) "collapse" of the statevector. 
This is the replacement of a statevector equal to a sum of vectors (each describing a different outcome of 
p% ' an "experiment") by one vector in the sum. One might guess that this dual evolution is indicative of a 

^ . fundamental deficiency in present day physics. In the hope of finding new physics, one may begin by trying 

to modify Schrodinger's equation so that the statevector undergoes only a smooth evolution, giving both the 
usual quantum behavior and the collapse behavior. 

This program, begun three decades ago,"'^'^ received a crucial impetus one decade ago from the work of 
Ghirardi, Rimini and Weber (GRW)^ and has evolved into the Continuous Spontaneous Localization (CSL) 
model. ^ At present, this is the only fully developed nonrelativistic collapse model, with definite predictions 
applicable to any nonrelativistic experimental situation. 

In the CSL model, a term which depends upon a randomly fluctuating field u;(x, t) is added to 
Schrodinger's equation. (The physical nature of this field is unspecified, but metric fluctuations,^ possi- 
bly with a tachyonic spectrum,^ have been suggested.) The probabilistic behavior of Nature is explained 
as due to our lack of control of the field w. When an experiment is under way,the particles in the sys- 
tem+apparatus interact with the particular sample field w that is present, causing a rapid evolution of the 
statevector to one of the alternative outcomes of the experiment. A different sample field leads to a different 
outcome. CSL also specifies the probability that a particular sample field w{x,t) actually occurs. When 
all possible fields are taken into account, together with their probabilities, the result is that each outcome 
occurs with (essentially — see next paragraph) the probability predicted by SQT. Thus two relations, the 
Modified Schrodinger Equation and the Probability Rule constitute the CSL model. 
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As with any modification of SQT, one expects and liopes for certain specially designed experiments 
where SQT and CSL lead to different predictions, making tests possible. For example, one such test, presently 
not practicable, is a two slit interference experiment with a sufficiently large bound state object. Once the 
two wavepackets for the object leave the slits, SQT says that their amplitudes will never change so that 
interference is possible at any time. CSL says that the amplitudes will fluctuate and, after a long enough 
wait, eventually one of the packets will become negligible in amplitude, giving no interference pattern. Such 
an interference experiment with e.g., 90A diameter drops of mercury^ over a time interval of seconds could 
provide such a test. However, this is a difficult experiment. For example, it is hard to prevent the two 
packets from being put into different angular momentum eigenstates by interaction with the environment, 
and then they would not interfere for this reason.*^ 

The tests that are most practicable at present stem from the consequence of CSL that the collapse 
process imparts energy to particles. (One may think of this energy as provided by the field w.) The reason is 
that collapse entails the narrowing of wavefunctions. By the uncertainty principle, this leads to an increased 
momentum spread, and thus to an increased energy. Thus, any bound ground state of, e.g., atoms or nucleii, 
will be excited by the collapse part of the Schrodinger equation. ^'^'^ The usual part of the Schrodinger 
equation will describe the radiation emitted as the system returns to the groimd state. Also, a free charged 
particle is "shaken" by the field w, and so it will radiate^". Similarly, the quarks inside a proton should 
be excited, and the proton should radiate mesons. Thus, a signature of CSL is that matter should emit 
"spontaneous radiation." In section 3, various radiation rates arc given. 

It is interesting that, at this time, quite a number of low noise experiments are being undertaken to look 
for radiation appearing in an apparatus for a variety of reasons, e.g., because of collisions with purported 
Dark Matter. Some of these experiments are sensitive enough to provide useful constraints on the parameters 
of CSL. 

Some experiments which look at radiation appearing in a slab of Germanium are described in section 
4. In section 5 we show how data from one such experiment, "Rico Grande," applied to the spontaneous 
excitation of a proton in a Ge nucleus, can provide a limit on the relative collapse rates of neutron and 
proton. Only a rough estimate is given because greater precision requires a more careful calculation of Ge 
nuclear dipole matrix elements than we are prepared to give here. 

In section 6 we consider the spontaneous ionization rate of a Is electron in a Ge atom. In a previous 
paper, it was argued that the upper limit on this rate given by a single data point from the "TWIN" 
experiments^ suggested that the coupling of the field w to an electron or nucleon (here assumed to be the 
same for neutron and proton) is proportional to the particle's mass, supporting previous proposals that there 
is a connection between gravity and collapse^'^^'^^. Unfortunately it subsequently turned out that the data 
point was inaccurate (see section 4). However, the more complete analysis on "COSME" data presented 
here gives essentially the previous result. 

Eventually, from such experiments, one may hope that the collapse rate parameter of CSL will either 
be constrained to be so small that the model will be ruled out — or that spontaneous radiation from collapse 
will actually be observed! 



2. CSL 

Underlying CSL are two mechanisms. One is the Gambler's Ruin mechanism. This explains how the 
random process embodied in the noise w produces the probabilities of SQT for the collapsed states. The 
other is the GRW "hitting" mechanism. It allows collapse to occur rapidly, for macroscopic objects, to states 
which we see around us (localized objects), while microscopic objects are scarcely affected.^ 

Here is the Gambler's Ruin analogy. Suppose, at the beginning of a game, gambler 1 (2) starts with 
<^i(0) ((^2(0)) dollars, and di(0) +<i2(0) = 100. This is to be analogous to the initial statevector \ip{0) >= 
ai(0)|l > +02(0)12 >, with the correspondence di{0)/WO |ai(0)p. The gamblers toss a coin (analogous to 
the fluctuating fleld w) and, depending on the result, one gives a dollar to the other. As the game proceeds, 
the di{t) fluctuate, just as do the squared amplitudes |ai(<)p. One gambler finally wins all the money and 
the game stops, with e.g., gambler 2 winning with probability (i2(0)/100. Precisely analogously, collapse 
finally occurs, e.g., with \tlj{t) >— > 0|1 > +1|2 > with probability |a2(0)|^. This is, of course, the probability 
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predicted by SQT of the outcome |2 > if |1 > and |2 > represent two states of an apparatus. 

The GRW model postulates a physical process which produces a sudden random change ("hit") of 
a many-particle wavefunction: the wavefunction is multiplied by a gaussian function exp— (x„ — z)^/2a^, 
where x„ is the position coordinate of the nth particle. The center of the gaussian, z, is chosen according 
to a Probability Rule which depends in a certain way upon the wavefunction, making it most likely that 
z is located where the wavefunction is largest. The effect of a hit is to narrow to width a the part of the 
wavefunction which depends upon the nth particle: GRW chose the mesoscopic length a « 10^"'' cm. A hit 
on one particle occurs at a slow rate A: GRW chose A ~ 10~^^sec~^ « once in 300 million years. But, each 
particle is equally likely to be hit, so a hit on an N particle object occurs rapidly, on average in 1/XNsec. The 
wavefunction of an N particle object in a superposition of states, each describing the object in a different 
place, has the particles entangled in such a way that one hit on one such particle causes the wavefunction 
to collapse in 1/ XNsec to one of the states in the superposition. This explains how macroscopic objects are 
always observed as localized. (A defect of this model is that the (anti) symmetry of the wavefunction is 
destroyed by the hitting process.) 

CSL may be thought of as embodying a continuous hitting process. A hit occurs every At sec, but the 
hit wavefunc;tion is multiplied by At and added to the original wavefunction. Thus the wavefunction shape 
fluctuates gradually (not suddenly as in GRW's model), and the gambler's ruin dynamics is obtained. The 
(anti) symmetry of the wavefunction is properly preserved in CSL. 

The CSL modified Schrodinger equation is 



d\ip,t >^ 
dt 



= -iH\^,t>y, j dx[w(x,t) -2A^(x)]2|^,t >^ (2.1) 



Given an arbitrary field 'w{x,t), one may solve (2.1) to find how the statevector evolves under its influence. 
Operator A(x) is 

(x - z)^ 

^(x) = V5„-4-T /rfziV„(z)e 2a2 (2.ia) 

The integral in (2.1a) is essentially the number of particles in a sphere of diameter a centered around x. 
Na{z) is the number density operator for particles of type a, e.g., electrons , protons and neutrons, and Xg"^ 
is their one-particle collapse rate (see Eq. (2.5) below). 

Thus, the parameters which characterize CSL are a, X and the ratios of the 5a s, e.g., for ordinary 
matter, gn/dp and ge/gp where the subscripts n, and e refer to the proton, neutron and electron. With 
no loss of generality we may take = 1, so A is an individual proton's collapse rate. 

CSL requires a second equation, giving the probability density P{w) functional that the field w(x, f) 
occurs: 

P{w)= ^<il!,t\tij,t>^ (2.2) 

Because the evolution (2.1) is nonunitary, the norm of the statevector >^ changes with time. Eq.(2.2) 
says that the most probable fields to occur are those which lead to statevectors of largest norm. 

Eqs.(2.1) and (2.2) ensure that an initial statevector, which is in a superposition of states of different 
particle number density, evolves toward one of these states for each probable w — and that the ensemble of all 
evolutions is such that each final state occurs with (essentially) the SQT probability. However, we shall not 
discuss here how Eqs.(2.1), (2.2) lead to collapse for an individual statevector. Instead we shall go right to 
the appropriate object for discussing experimental predictions, the density matrix p. The evolution equation 
for the density matrix whch follows from Eqs.(2.1), (2.2) can be shown to be * 

d < x\p(t)\x' > . ..^^ , ,,. , 
= -i< x\[H,p{t)]\x' > 

^ Af w (2.3) 

~ 2^212 SaU)9a{k)[^{^3 - Xfe) + $(x; - x;,) - 2$(x, - X^)] < x\p{t)\x' > 

j=l fe=l 
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where \x >= |xi,X2, . . . > is the position eigenstate for all particles and 

$(z) = e 4a2 (2.4) 

As a simple example of how Eq. (2.3) works, set H = so as to concentrate on the collapse dynamics 
alone, and consider a clump of particles of type a in a superposed state. That is, let the initial state be 
ai(0)|l > +a2(0)|l >, where |1 > and |2 > each describe N particles of type a in a localized state with 
dimensions << a, but with centers of mass of the two states at a distance >> a apart. Then $(xj — Xfc) w 1 
if Xj, x;j are both located in region 1 (or both in 2) and $(xj — xj.) « if Xj, x^ are located in regions 1 
and 2 respectively. Therefore, Eq.(2.3) yields 

^^''gf = E E(l + 1 - 2 • 0) < l\p{t)\2 >= -XglN^ < l\p{t)\2 > (2.5) 

J = l k=l 

showing that the off-diagonal elements of p decay at the rate AA'^^. This illustrates how the collapse rate is 
large for a superposition of states describing a large number of particles in different locations. 



3. Excitation Rate Predictions 

As mentioned in section 1, a byproduct of the collapse process is that particles gain energy. It is easy 
to show, using Eq.(2.3), that the average total energy H{t) —"^^^^ '^^ 
increases according to 



TrHp{t) {H = Er=i Pi/2Mi + ^(xi . . .xa,)) 



dH{t) _ 3A 
dt 



N 



(3.1) 



Assuming = 1 for all particles, and using the GRW values for A and a, then 10-^'* nucleons gain « .3 
eV/sec and 10^^ electrons gain « 600 eV/sec. This is quite small, corresponding to a temperature increase 
over the age of the universe of « .GG1°K and 2°K respectively. (The low particle density in the universe 
assures that the effect of this increased energy on the the cosmic radiation bath is negligible). If we take 
10^* electrons as roughly the number in a cc. of condensed matter, this corresponds to an energy increase of 
about lO^-'^"'' joules/scc, which is close to the experimentally detectable lower bound by present day bolometric 
measurements.^^ This is much less sensitive than the experiment discussed here. 

However, while (3.1) gives the average behavior, there are infrequent but large energy fluctuations. 

The energy increase in Eq.(3.1) is the sum of increased internal energy and of increased center of mass 
energy ifem = ^IJ^M + V,m{Q) (Pcm = Ef=i Pj, M = Ef=i M,-, Q = Ef=i M,x./M). From Eq.(2.3) 
we find 



dt 



.2 N N 

^ EE 

j=i k=i 



dHcm{t) 3A?l^ \ ^ \ - _1 



(Xj - Xfc) 



_3A^ 

4a2M 



N 



E^«o) 

■i=i 



6a2 



(3.1a, 6) 



where the first term in the expansion (3.1b) dominates if the system under consideration, like an atom or 
nucleus, has dimensions « a. The condition for the total energy increase to be completely due to the center 
of mass energy increase to order a~^, i.e., for there to be no internal excitation to this order, is found by 
equating (3.1) to (3.1b): 



N N -.2 N N y 
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i.e., if ga{j) = CMj. Therefore, for such mass proportionahty of the couphng constants, the internal energy 
does not increase to order a~^. Moreover, the leading term in the internal energy increase, proportional to 
a"^, is compensated by an identical decrease in the center of mass energy since, by (3.1), the total energy 
increase vanishes to order a"^ and higher. 

If the coupling constants are not mass-proportional, the internal excitation rate ~ a"^ is found as 
follows. Consider a transition from a state \ip > \x > to a, state > \x' >; where \tp > is an initial bound 
state, \(f> > is an orthogonal final state, \x > is an initial state of the center of mass and \x' > is an arbitrary 
final state of the center of mass. The probability per second of a transition from \tp > to \4i >, regardless of 
the final center of mass state is P = S|x'> ^ ^'1 ^ ^IpWI^' > \x' >> where p(0) = |x > l^^" >< V'l < xl- 
Expansion of Eq. (2.3) to first order in (dimension of system/a)^ yields^i 

A = 2^ < > ■ < m\<t> > (3.2) 

where R = Ylf=i 9aU)'^j, = - Q. 

It is the predictions of Eq.(3.2) that we shall test in this paper. It is worth remarking that the matrix 
element in Eq.(3.2) involving a charged particle's is the same as that involved in describing an electric 
dipole transition between lip > and >. Thus one can evaluate this contribution to (3.2) either by 
calculating the relevant matrix element or by expressing it in terms of measurable transition rates. Indeed, 
as we shall see in sections 5 and 6, using R = if 5^ ~ M^, it is possible to express the matrix elements of 
one type of particles in terms of another type so one need only calculate or measure the matrix elements of 
the excited particle type to evaluate (3.2). 

Incidentally, by summing Eq.(3.2) over all states \(p > orthogonal to \ip >, we obtain P^, the total 
probability/sec for excitation of |V' >: 

A^ = ^<^|[R-<^|R|'A>]'I^> (3.3) 

Assume the GRW values for A and a. For an atomic electron undergoing spontaneous excitation from, 
e.g., the Is state of an atom with atomic number Z to a higher energy state, bound or free, the order of 
magnitude of A is « gplO^^'^/Z^scc^i. For, e.g., a proton in an outer shell of a nucleus of mass number A 
it is w g^lO^'^^A^^^sec^^. With such rates, the Is electrons in lO-^'* such atoms would be expected to provide 
Ri g'^lO/Z'^ photon pulses each second while each spontaneously excited proton in 10^^ such nucleii would 
be expected to provide w g^.3A^^^ gammas each year. This large difference in rates explains why we are 
able to obtain good experimental limits on the electron's coupling constant ge in section 6, but not on the 
neutron's coupling constant gn in section 5. 

Although we shall only apply Eq. (3.2) in our data analysis, for completeness we include Eqs. (3.4), 
(3.5) below which could be appled if more accurate data becomes available. We note cigain that if is 
mass-proportional (i.e., if for protons gp = 1, then for neutrons Qn ^ 1.001 and for electrons Qe ^ .00054), 
then R = 0, and therefore (3.2) vanishes. We should then need P to order a~^: 

3 3 

= 1^1' < "^'^l^ > I' + 2 E E I < <^l'5""IV' > I'] (3-4) 

m=l n=l 

where 5™" = J2jLi 5a(i)(R'i)'"(R'i)") = J2n=i corresponds to a charged particle, its matrix 

elements here describe electric monopole or quadrupole transitions. These matrix elements are smaller than 
the corresponding matrix elements in Eq.(3.2) by the factor (size of bound state/a)^. 

Fu^^ has considered the spontaneous electromagnetic radiation of a free charged particle in CSL, ob- 
taining the probability /sec/energy of radiating a photon of energy E = %k: 



dE ^"i^^l^J E " -eT^ counts/sec/keV (3.5a, 5) 

(the infrared divergence is treated as usual). Eq.(3.5b) gives the rate for an electron with ge = 1. This 
radiation rate (3.5) for free particles is smaller than the excitation rate (3.2) for bound particles by the 
factor e^/hc w 1/137 and by the replacement of (bound state size/a)^ by (Compton wavelength/a)^. We 
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note, with mass-proportionality, that the spontaneous radiation rate for free electrons is the same as for free 
protons, on account of the factor (jjal^a)^ in (3.5a). 

This completes our collection of equations giving CSL spontaneous excitation rates. 

Wc now consider the spontaneous excitation of valence nuclcons in a Germanium nucleus and the 
spontaneous ionization of Is electrons in a Germanium atom, using data from what are, at present, the 
lowest noise relevant experiments. 



4. Experiments. 

The data used in the analysis-'^* of the atomic excitation comes from a small (253 g) p-typc coaxial HPGe 
crystal, "COSME", built specifically for a Dark Matter search^^. It features, at only 1.6 keV, the (so-far) 
lowest energy threshold of any detector dedicated to such searches, and has as well an excellent resolution 
of 0.43 keV (FWHM) at 10.3 kcV. Special low-radioactivity measures were taken, such as mounting the 
detector on a specially-designed electroformed copper cryostat, shielding of electronic components close to 
the crystal with 450-yr-old lead, and use of a 2000-yr-old roman lead layer in the innermost part of the 
shielding. Additional photon, neutron and vibrational shielding were used. 

The detector set-up was installed in the Canfranc-1 underground laboratory in the Spanish Pyrenees, at 
a depth of 675 meters of water equivalent. The microphonic component characteristic of very low-threshold 
detectors, extending up to « 15kcV, was filtcrcd-out using specially-developed techniques While the low- 
energy background level was slightly higher than that in the "TWIN' detectors the improved resolution 
— typically inversely proportional to the mass of the crystal — allows one to impose more stringent limits on 
a sharply defined signal that might be buried in an otherwise featureless background, as is the case for the 
emitted radiation in CSL. 

At this point a remark is in order: unfortunately the TWIN data used to extract CSL limits on in 

reference 12 belonged to a preliminary set coming from un-amplificd digitized pulses. Later comparison with 
the spectrum collected with a multichannel analyzer showed that this earlier data was corrupted at energies 
below « 200 keV, i.e., a large fraction of events were not recorded. This faulty set was not used in other 
TWIN results, namely for double-beta decay ^ Dark Matter "^^ or electron half-life In reference 12, 
only one data point was used, the (erroneous) rate .049 counts/keV/kg/day at 11 keV (where the highest 
predicted value of spontaneous radiation occurs) to set an upper limit on g^. In the present work, it is 
necessary to perform a Chi^ analysis of a fit to the whole predicted shape of the spontaneous radiation, in 
order to obtain a limit comparable to the (erroneous) limit of reference 12. 

The data used in the analysis of the nuclear excitation comes from the "Rico Grande" crystals, which are 
part of the IGEX ensemble of large enriched germanium detectors^'^, dedicated to searching for neutrinoless 
double-beta decay. The two detectors from which these data were extracted have a fiducial mass of «2 kg 
each and are enriched to 86% Gerg and 14% Ge74. As a result, the prevailing source of background in the 
energy neighborhood of interest here is this two-neutrino double-beta decay from Geyg which, however, cuts 
off at an energy below the region employed in the analysis in section 5. At the time of collection of the 
present data, the Ricos were operated in the Homestake mine in similar conditions to TWIN. 



5. Constraint on gn/ gp^ 

We now apply Eq. (3.2) to the spontaneous excitation rate of a single proton or neutron in a Ge nucleus. 

The Jparity for the ground state of Ge is 0+ for the even-even nuclides Geyo (20.6%), Ge72 ( 27.4%), 
Gc74 (36.7%) and Ge76 (7.7%)), and it is 9/2-1- for Ge73 (7.7%). Spontaneous excitation is predicted (the 
matrix elements are nonvanishing) for a transition from the ground state to 1- states in the case of the 
even-even nuclides and to 11/2-, 9/2- or 7/2- states in Ge73. Prom the point of view of the shell model, 
among other possibilities, a proton or neutron can be excited from its ground state valence level to a higher 
energy state. 

The return of a proton from the excited state to ground could be direct, via an electric dipole transition. 
It could also proceed indirectly, through intermediate states or internal conversion and, in the case of a 
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neutron it must proceed indirectly, through magnetic transitions. In any case, the lifetimes are in most cases 
so short that a photon pulse would rapidly appear at the energy difference of the two states. Therefore, by 
looking at the data for a signature peak of instrumental resolution width at the expected energy, one may 
hope to observe the radiation resulting from these spontaneous transitions or at least get an upper limit on 
their rate. 

Not only does the matrix element of Rj for the excited particle not vanish, but the matrix element of 

Rj for the other particles also does not vanish due to their dependence on the center of mass operator. We 
can relate the matrix element of the protons to that of the neutrons by using S^Mj-Rj = 0, which implies 
that 



E^-^R™ = -(Mp/M„)Sf^iRp, (5.1a) 
(we neglect the electron contribution of o(Me/(M„)) so, setting Mp/M„ = 1, the matrix element in (3.2) is 

< (A|5pS,llRpi + gn^tJi^Rnil^ >= [1 - Qn] < <^|S,iiRpi|V > (5.16) 

(remembering t/p = 1). Thus, from Eqs. (3.2) and (5.1), the excitation rate F in sec~^ of e.g., one of the 
four valence protons from the ground state \tp > to one of the 12 degenerate excited states \<j) > (any of the 
four protons may be excited, and the 1- state of this proton plus its subshell partner can have three possible 
orientations) can be expressed purely in terms of proton matrix elements: 

r = A[i _ g„]2i2| < <^|Sf^,Rp# > p 

I (5.2a, 6) 

=^[1 -5n]'i2S„^| < <A|s£fi?piyi,„(^,,0,)|v > p 

A similar expression may be written for the excitation rate of a neutron totally in terms of neutron matrix 
elements. 

It is worth remarking that the expression for the lifetime r(0 ip) of one of the excited state |(/) > to 
decay by an electric dipole transition to the ground state \ip > can be written in terms of the same matrix 
elements as appear in (5.2)^^: 



1 _ 16ttc / E \3 f ' 
T{(t) V) ~ 9 

where E is the energy difference of the two states. Thus we may express F in terms of this lifetime: 



' 3 2 



A . ,2 9 /?ic\3 1 
2^^ ^^"J c(e2/M \e) t(0 ^ V) ^ ^ ^ 

Eq. (5.4) would be useful if we have the experimental lifetimes and branching ratios of the state </). Unfor- 
tunately, in this case we do not, so we are forced to estimate the matrix element in (5.2b). 

In this paper we shall approximate the matrix element by using the same "very rough estimate" employed 
by Blatt and Weisskopf^^ for calculating the lifetime r. They assume the radial wavefunction of the proton 
in both states (j) and V is 6(J?o - R)[Z/RoY^'^ where 6 is the step function and Ro = 1.4 x 10~^^A^/^ is 
the nuclear radius. As they point out, the actual radial integral is expected to be "somewhat smaller" since 
radial wavefunctions oscillate: say, /? times smaller. We obtain from (5.2b) the result 

F = (5 5) 

(a numerical factor 9/8 has been replace by 1). 

The expected total count C for an experimental run of D kg-days from a transition due to a nuclide 
which comprises the fraction X of the 8.3 x lO-^'* atoms/kg in common Ge is found from (5.5), with the GRW 
parameters, to be 

C=^[l-gr,fXD (5.6) 
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For example, consider a transition in Ge74 to the 2165 keV 1- state. We shall use the data from the Rico 
Grande experiment^^^ similar to COSME but with D = 1.135 kg-yrs=414.3 kg-days (COSME's D = 85.24 
kg-days) with X = .14 (COSME's X = .37). Denoting by 

Cexpt the upper limit on the number of observed 

counts, we obtain 

1 1 
1 + .26^^ > 5n > 1 - .26-^ (5.7) 

For this experiment, the upper limit (obtained from the counts under a Chi^ fit to a background quadratic plus 

the expected experimental resolution shape centered on 2165 keV) is Cexpt = -89 counts at the 68% confidence 
level (3.9 counts at the 95% confidence level). For /? w .3 one obtains 1.8 > gn > -2. However, our choice 
of /3 is just hypothetical, as we have not made the effort to seriously evaluate the matrix element^^'^''. The 
points to be made are that the range of g„ is not so far from g-p = 1 and that the various numbers involved in 
calculating (5.7) tantalizingly contrive to be on the edge of showing mass proportionality. Indeed, this would 
more easily be shown with a larger value of A/a^ than the GRW value: for instance, with X/a^ = 100 X/ a Qj^y^r, 
the above inequality becomes 1.1 > ,g„ > .9. But, with the GRW parameters, it would require a long counting 
time and a proper calculation of matrix elements before one might say that gnl 9p ~ 1- 



6. Constraint on ge/gp- 

Here we apply Eq. (3.2) to calculate the spontaneous ionization rate of the Is electrons in a Ge atom. 

If a Is electron is spontaneously ionized, the remaining electrons in the atom rapidly cascade into the 
(singly ionized) ground state, emitting a photon pulse of 11.1 keV (the ionization energy of a Is electron). 
The ionized electron also deposits its kinetic energy in the Ge sample, which augments the energy of the 
pulse. Thus the signature of these events is a distribution of photon pulses of energy E > 11.1 keV. 

In reference 12, a Hartree calculation of the matrix element in (3.2) for the electrons was numerically 
performed, where \ip > is the ground state of Ge and |0 > is a state where a Is electron is ionized. The result 
of the calculation may be expressed as a function C{E) which gives the expected pulse counting rate if GRW 
parameters are assumed and if the electron is totally responsible for collapse (i.e., ge = 1, gn — gp = 0). 
C{E) is zero for E < 11.1 keV, abruptly rises to 5370 counts/keV/kg/day at -E = 11.1 keV, and decays 
in roughly exponential fashion, with the value « 4000 counts/keV/kg/day at E = 12 keV, and « 1500 
counts/keV/kg/day at S = 16 keV. 

If we assume that g„ = = 1 (as we shall hereafter do) then, as in the preceding section's Eq. (5.1), 
we can express the matrix element for the nucleons as —Me/Mp times the matrix element for the electrons. 
Putting this into Eq. (3.2), the resulting rate T in counts/sec/kg/day may then be written as 



(A/a2)Gw Mp 

We note that the rate in (6.1) vanishes if there is mass proportionality {ge/gp = M^/Mp). 

Figure 1 shows a graph of the counts/. lkeV/85.24kg-day from COSME in the energy range 5 to 17 keV. 
A recent paper^^ describes a search with the same apparatus in a similar energy range. These authors were 
looking at the TWIN data for a signature 11.1 keV peak resulting from a hypothesized violation of charge 
conservation (in which a K-shell electron decays to neutrals, and the other electrons in the atom readjust). 
They fit the data to three x-ray peaks (Cu, Zn and Ga) known to result from cosmogenic excitation of the 
Ge isotopes in the sample under observation, together with a background quadratic polynomial plus the 
hypothesized process, and they look at the difference between the fit and the data at the 68% and 90% 
confidence levels. 

We employ here the same procedure. Fig. 1 (solid line) shows the best fit to the data by this method, 
without the hypothesized process, a multiple of C{E). Superimposed upon this fit is a graph of lO^'^C(-E) 
(dash-dotted line), folded in with the detector resolution shape (a gaussian of standard deviation .18 keV). It 
is clear that this is by no means a good fit to the data, so the coefficient of C{E) in Eq. (6.1) is considerably 
smaller than 10~^. 
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Also shown in Fig. 1 is the fit with the hypothesized process at 4.2 x 10~^C{E) (dotted line) which 
corresponds to the 68% confidence leveL Not shown is a similar shaped curve at 9.7 x 1Q~^C{E) which 
corresponds to the 95% confidence level. 

If we assume that the parameters A and a are the same as those given by GRW, then we conclude from 
Eq.(6.1) that, at the 68% confidence level, [c/e/ffp - Me/Mp]^C{E) < 4.2 x 10-^C{E) or 

< % < 13^ (6.2) 

- 9p- Mp 

Thus, according to CSL with the GRW parameters, nucleons are mostly responsible for collapse. 

It is worth noting that it would take an improvement in the experimental limit by e.g., a factor of 1/300, 
which results in .SM^/Mp < gel dp < l-7Me/Mp, to suggest that ge/dp = may be ruled out. However, it 
should also be noted that we need not be wedded to the GRW parameters. Thus an increase of A/a^ by a 
factor of 300 or more would have the same effect. On the other hand it would take a decrease in A/a^ by 
a factor of 4 x 10~^ or more for the limit obtained in this experiment not to suggest that nucleons collapse 
more rapidly than electrons. 
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Figure Captions 

Figure. 1 A graph of the COSME data in the region 5-17 keV is shown along with the best fit to the 
three known X-ray peaks plus a quadratic polynomial background (solid c;tu-vc). Two additional curves are 
shown, corresponding to predicted rates folded in with the experimental resolution (a gaussian of width 
.18keV). The dash-dotted curve corresponds to 10~^C{E) and the dotted curve to the 68% confidence level 
value of 4.2 x 10-^C{E). 
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